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Abstract 

We derive radial and vertical distributions of HI and H2 gas densities in our 
Galaxy by using the terminal velocity method. We calculate the molecular frac- 
tion (/moi) defined as the ratio of the molecular hydrogen to total hydrogen gas 
density at galactic longitude I = 33° ~ 64° and galactic latitude b = —2° ~ +2°. 
The thickness of the molecular dominant region (/moi > 0.8) is approximately 
constant (109 ± 12 pc) at galactocentric distance R ~4.7— 7.2 kpc. The molecu- 
lar fraction decreases suddenly at a critical height from the galactic plane, below 
which the gas disk is almost totally molecular, while it is almost atomic beyond 
this height. We show that the vertical /moi variation can be reproduced by a 
model which takes into account the phase transition between HI and H2 gases 
in the interstellar matter. 

Key Words: Atomic Hydrogen — the Galaxy — Interstellar Matter — Molec- 
ular Hydrogen 
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1 Introduction 



The major constituents of the interstellar matter (ISM) in the Galaxy are the neutral 
atomic (HI) and molecular (H2) hydrogen gases. The distributions of the HI and H2 gases 
have been obtained through observations of the 21-cm and CO-hne emissions, respectively. 
It is known that the HI gas is distributed in the outer region of the Galaxy, while H2 gas is 
more concentrated toward the galactic center (e.g.. Burton 1988; Combes 1991). However, 
the mutual relation of the distributions of HI and H2 gases in the Galaxy, which will be 
deeply coupled with the transition between these two phases of the ISM, has been not yet 
investigated in detail. 

Elmegreen (1993) has proposed a phase-transition theory of the HI and H2 gases, in 
which the molecular fraction (/moi)j as defined by the ratio of the H2 gas mass density to 
the total hydrogen gas density, is determined by the interstellar pressure, UV radiation 
intensity and the metal abundance. Sofue et al. (1995) and Honma et al. (1995) have 
studied the radial variation of /moi in several galaxies, and found that the molecular 
fraction decreases suddenly at a certain galactocentric distance, which they called the 
molecular front. They further modeled the molecular front phenomenon based on the 
phase-transition theory of the HI and H2 gases in the ISM. 

Since the Milky Way Galaxy has been observed extensively both in the HI and GO 
lines at much higher linear resolution than external galaxies, we may be able to study the 
molecular front phenomenon in the Galaxy not only in the radial direction but also in the 
perpendicular direction to the galactic plane. In this paper we analyze latitude-velocity 
(b — v) diagrams of the HI- and CO-line emissions of the Galaxy by applying the terminal 
velocity method, and derive the HI and H2 gas distributions in the {R, z) plane, where 
R is the galactocentric distance and z is the height from the galactic plane. We further 
derive the distribution of the molecular fraction in the (i?, z) plane, and investigate the 
hydrogen phase transition in the z direction. We also model the result based on the 
phase-transition theory of ISM. 

2 Data and Procedure 
2.1 Data 

We use the HI survey of Bania and Lockman (1984) observed with the Arecibo 304-m 
radio telescope. The HI data are in the forming h — v diagrams, and maps are sampled 
at intervals of Afe = 2' with a half-power beamwidth of 4', and the latitude coverage is 
±3°. We also use the GO survey of Knapp et al. (1985) observed with the 7-m telescope 
at the Bell-Telephone Laboratories. The GO maps are sampled at intervals of A6 = 2' 
with a half-power beamwidth of 100", and the latitude coverage is ~ ±2°. The longitude 
coverage of the HI and GO data are Z = 33° ~ 64° and Z = 10° ~ 77°, respectively. Both 
maps are spaced at equal intervals in longitude: AsinZ = 0.025. We used the Arecibo 
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HI data, which have the highest resolution for HI observations, in order to compare with 
the CO data of a resolution 100" and to examine thevertical (latitude) variation of gas 
density. Table 1 shows the observational parameters of the data, and Figure 1 shows 
examples of the HI and CO b — v diagrams superposed on each other. 

- Table 1, Figure 1 - 



2.2 Molecular fraction 

We adopt the following relations between the column density N of the atomic and molec- 
ular hydrogen and the intensity / of the HI and CO line emissions: 

iVn [cm-2] = Ci/hi (1) 

and 

7VH,[cm-2]= C2/C0, (2) 

where Ci and C2 are the conversion factors. In this paper we regard C as a constant and 
adopt 

Ci = 1.8 X 10^^ cm-V(Kkms-^) (3) 

and 

C2 = 2.8 X 10^° cm-V(Kkms-^) (4) 

(Burton 1988; Bloemen et al. 1985). However, Arimoto et al. (1996) point out that the 
CO-to-H2 conversion factor in spiral galaxies is dependent on the radial distance from the 
Galactic Centeras 

C2 = C°10°=^^^/^% (5) 

where the coefficients have been given as = 0.9 x 10^° cm~^/(Kkms^^) and Re = 6.6 
kpc. Although our study will be based on an assumption of a constant C2 factor as 
equation (4), we examine the effect of the radial variation of C2 on our result, and discuss 
it in the last section. 

The number density n is given by 

n = N/Ar, (6) 

where Ar is the depth along the line of sight. 

We define the molecular fraction /moi as the ratio of the molecular hydrogen gas density 
to the total hydrogen gas density, which is expressed by mass density p or number density 
n by 

n _ Phz _ 2nH2 

/mol , I o ' \ ) 

PR + Ph2 nu + 2 nii2 
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2.3 Tangent points and Terminal-velocity method 

Assuming a circularly symmetric model of the Galaxy, we can obtain distances to the 
emission at terminal velocities in the first and fourth quadrants (0° < / < 90° and 
270° < / < 360°). The radial velocity at each galactic longitude / in the first quadrant is 
given by 



Vr — 



^■V{R) - Vo 



sin I, (8) 



where V{R) is the rotation velocity at galactocentric distance R, Ro{— 8.5 kpc) is the 
distance of the sun from the galactic center, and Vo{— 220 kms~^) is the rotation velocity 
at J? = Rq. 

Terminal- velocity emission comes from the tangent point at a galactocentric radius 

R^Ro\sml\, (9) 

and the galactic latitude b corresponds to the height from the galactic plane as 

z = Rq cos / tan 6. (10) 

The actual velocity profiles do not have a sharp cutoff due to the velocity dispersion of 
the gas. So, we regard the integrated emission within a certain velocity range of width Av 
around the terminal velocity as a contribution of the gas at the tangent point. We adopt 
a velocity width Av = 6.5 km s^^ at any longitude, which we found to be a typical width 
of the CO velocity profile corresponding to a terminal-velocity emission component at 
6 ~ 0. In fact, this value is consistent with the velocity dispersion obtained by Malhotra 
(1994), who finds that the value of the velocity dispersion varies 3.8 km s^^ at 2.5 kpc to 
7.1 kms~^ at 7.5 kpc. The depth along the line of sight corresponding to this velocity 
width is given by 

1/2 

(11) 



Ar = 2i?o sin I 

0.50 i?o sinZ. 



Vq-Av 
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Based on these relations we calculated number densities of H and H2 at each longi- 
tude and latitude using latitude-velocity {b — v) diagrams, and obtained vertical density 
distributions as a function of height z from the galactic plane. Then, we arranged the 
thus obtained z distributions at various longitudes in the order of the distance from the 
galactic center, and obtained two dimensional maps of the distributions of H and H2 gases 
in the {R, z) plane. 



3 Result 

3.1 HI and H2 Distributions 

As seen from the b — v diagram (Figure 1), the HI gas is distributed smoothly in the 
velocity direction at any longitude, and therefore, at any point along the line of sight. 
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On the other hand, the H2 gas distribution is clumpy. The HI and H2 distributions in 
the z direction are also clearly different (Figure 2): the HI distribution has a ffat plateau 
near the galactic plane and gentle skirts, while most of H2 profiles have multiple peaks, 
and the maximum density is much higher than that of HI. The CO distribution is not 
centered on the Galactic plane, which is consistent with previous studies (Malhotra 1994; 
Dame et al. 1987). 

- Figure 2 - 

Figures 3a and h show the {R, z) distributions of HI and H2, respectively. These figures, 
therefore, represent cross sections of our Galaxy along the tangent velocity circle. The 
density distribution of HI gas in the z direction extends for a few hundred parsecs, while 
the extent of H2 is not more than 100 pc: there is little CO emission at high latitudes. 
Note, however, that these two cross sections are not fully sampled in the longitude direc- 
tion, because the observations have been obtained every A sin / = 1/40 (AZ = l.°4 / cosZ : 
Ai? = 212.5 pc) with a beam width 4' for HI and 100" for CO. 

- Figure 3 - 

In Figure 4, we show an (i?, z) distribution of H2 gas for the whole range covered by 
the present CO data from the galactic center region (/ = 10°, i? = 1.49 kpc) to the solar 
neighborhood (/ = 77°, i? = 8.29 kpc). This figure includes the region shown in Figure 
36. Figure 5(a) shows radial distributions of HI and H2 gases averaged at \z\ < lOOpc. 
A radial distribution of CO emissivity in our Galaxy has been obtained by Robinson et 
al. (1988) and Sanders et al. (1984). Our result is approximately agree with theirs. 
The variation of the H2 number density is much larger than that of HI. Note that the 
highest density region for H2 in Figure 4 and 5 corresponds to the '4- kpc molecular ring', 
and that the H2 gas density in the inner region within 3 kpc and outside of 7 kpc is an 
order of magnitude smaller than that in the 4-kpc ring. Figure 5(6) shows the same as 
Figure 5(a) but using the C2 factor as given by equation (5), which takesinto account the 
radial variation of the CO-to-H2 conversion factor. The H2 density distribution changes 
drastically in the inner region of the Galaxy, indicating much less molecular gas within 
Rr^ A kpc. 

- Figure 4, 5 - 

3.2 Variation of /moi 

We calculated /moi for a region shown in Figure 3, where both the HI and H2 densities are 
obtained. Negative values in the observed CO intensity (H2 density) have been replaced 
with zero. Figure 6 shows an example of the variation of /moi at I — 53°. It is evident 
that the value of /moi is almost unity near the galactic plane, and decreases rapidly to 
zero in a narrow range of 2; as 2; increases. 
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This phenomenon is commonly seen at all longitudes except at / = 61° and 64° (R < 7.5 
kpc) where the amount of H2 is small. This fact indicates the existence of a 2;-directional 
boundary of a region with a high H2 fraction, which we call the molecular front in the 
2;-direction. Most of hydrogen gas is in the molecular phase in the region near the galactic 
plane, sandwiched between the two molecular fronts, whereas in the outer (higher) region 
it is in the atomic phase. 

- Figure 6 - 

Figure 7 shows the two dimensional distribution of f^oi as a contour map in the {R, z) 
plane. The centroid of the molecular dominant disk, where /moi > 0.8, is waving, cor- 
responding to the fact that the molecular gas disk is corrugating (e.g. Malhotra 1994). 
It is remarkable that the thickness of the molecular dominant disk is almost constant at 
109 ± 12 pc. On the other hand, it is known that the scale height of the HI layer is larger 
than it. A thickness of HI layer between half-density points has approximately constant 
value: 220 kpc over 4 < i? < 8 kpc (Dickey and Lockman 1990). When we compare the 
thickness of molecular dominant disk with the scale height of HI gas, we can express that 
the atomic hydrogen layer wraps the molecular disk. 

- Figure 7 - 

Figure 8(a) shows radial variations of the molecular fraction calculated for number 
densities of HI and H2 averaged within \z\ < 50 and 200 pc of the galactic plane. In 
both cases the molecular fraction decreases from 0.8 to 0.6 with increasing R. Though we 
can not clearly see the phase transition in this region, curves are consistent with the fact 
that the molecular fraction is almost unity in the the central few kpc region (Sofue et al. 
1994). Particularly, the curve for \z\ < 200 pc indicates a decrease of /moi beyond R ~ 
7.5 kpc. This may be consistent with the existence of the radial molecular front indicated 
from an analysis of longitude- velocity diagrams of the Galaxy (Sofue et al. 1994). Since 
the H2 gas density over R > 8 kpc is much smaller than that of inner region (Dame et 
al. 1987), it is expected that the molecular fraction will decrease suddenly beyond this 
radius. Figure 8{b) shows the same variation but considering the radial dependency of 
CO-to-H2 conversion factor. The radial molecular front slightly change. However, no 
significant effect is found. 

- Figure 8 - 

3.3 Model Analysis 

The three crucial parameters which determine /moi are the interstellar pressure (P), the 
intensity of dissociative UV radiation field (C/), and the metal abundance (Z) which is 
related to the amount of interstellar dust and shields UV photon (Elmegreen 1993). Based 
on the phase-transition theory (Elmegreen 1993), Honma et al. (1995) has developed a 
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computing code to model the molecular front in galactic scale. By applying this code, we 
calculate the variation of molecular fraction in the 2;-direction. 

Since it is reasonable to assume that the pressure is proportional to the gas density, 
we consider the following two cases: first, P{z) is traced by the observed total gas density 
distribution in the z-direction (ntot), and second, P{z) is expressed by a gaussian function 
of z, considering the hydrostatic balance of the gas pressure with the gravity in the z- 
direction. According to van der Kruit and Searle (1982) we assume that the intensity of 
the radiation field is expressed as 

U{R, z) = Uq e^-^/^'^^sech^z/zo), (12) 

where Rjj — 5.0 kpc is a scale distance from the galactic center, and zq — 0.7 kpc is the 
scale height in the 2;-direction. For the metallicity, we do not have any information about 
the gradient in the z-direction. So, we assume that the value of Z is constant against z, 
while it varies with R as Z ^ Zq e-'^^-^'>/^^ . 

- Figure 9 - 

Figure 9 shows a thus calculated variation of the molecular fraction for / = 53°. The 
solid line corresponds to the observed ntot{z), and the dashed line to the gaussian distri- 
bution of P{z) in which the parameter of the gaussian function was taken so that it fits 
the observed ntot(^)- Even for such a round pressure distribution as a gaussian, /moi shows 
a shoulder-like variation, indicating the phase transition clearly. The curve of /moi calcu- 
lated for the observed ntot(^) shows a sharper variation at a critical height, and agrees 
well with the observed /moi variation (Figure 6) and reproduces the molecular front in 
the ^-direction. Though it has been considered that the most important parameter for 
the molecular front is the metal gradient (Honma et al. 1995), our results shows that the 
molecular front is reproduced without a metal gradient. 

4 Discussion 

The radial molecular front has been observed in our Galaxy as well as in many external 
galaxies (Sofue et al. 1994; Honma et al. 1995), in which the front is found at /? ~ 

5 — 10 kpc. In our result for our Galaxy, although /moi is almost unity in 4.8 < R < 7.2 
kpc, the radial decrease of /moi may occur at i? > 7.5 kpc, where the value of /moi drops 
from ~ 0.8 to less than 0.6. Since the HI gas density varies more gradually than the H2 
gas at R > 7 kpc, we may expect that the molecular gas is not any more dominant in 
the outer region of i? = 7.5 kpc. This indicates that the molecular front in our Galaxy 
is located at i? ~ 7.5 kpc and the transition region from H2 to HI and vice versa is as 
narrow as Ai? ~ 2 kpc in the radial direction. 

From the ^-directional variation of /moi, we found the molecular fronts in the direction 
perpendicular to the galactic plane. The molecular hydrogen disk is sandwiched by a 
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thicker HI envelope. In the molecular region near the galactic plane, between the two 
fronts, the ISM is almost totally in the molecular phase with /moi > 0.8, while in the 
outer region, the gas is almost totally HI. We stress that the transition between HI and 
H2 is taking place within the narrow range at a height of z ~ ±50 pc at < 7 kpc. 

In this paper, we have not taken to account a metallicity gradient in the z-direction, 
since the existence of the z-gradient of the metaUicity has not been reported. Even if it 
exists, the scale height of the variation would be a few hundred pc, because it will reflect 
the scale height of old population stars. Therefore, it will be reasonable to assume that 
the metallicity gradient in ^-direction is small or not present in our region of < 100 
pc. 

It has been found that the CO-to-H2 conversion factor (C2 factor) increases exponen- 
tially with the radius due to the metallicity gradient (Arimoto et al. 1996). This implies 
that the C2 value in the region considered here (from R=4 to 7 kpc) varies by a factor 
of 1.5. We have examined the radial variations of the H2 density and /moi by taking into 
account this effect. We found that the radial molecular front slightly moves outwards, 
but not significant change of the property of the /moi distribution in the ^-direction was 
found. 

We thank M. Honma for his valuable comments and for providing us with the code to 
model the molecular front. We are also indebted to Drs. T. M. Bania, G. R. Knapp, and 
S. Malhotra for providing us with the HI and CO data in a machine readable formats. 
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Figure captions 

Figure 1. A latitude- velocity (b — v) contour map of the HI brightness temperature at I = 39°. The 
ordinate is galactic latitude in degrees, and the abscissa is velocity relative to the LSR in kms~^. 
Contours are drawn at 10% interval of the peak brightness (114 K). Ab — v map of the CO emission 
at I = 39° is superposed in gray scale. 

Figure 2. Vertical distributions of the atomic hydrogen gas (doted line) and the molecular hydrogen 
gas (solid line) at I = 53°. The ordinate is the number density in HI cm~^ and H2 cm~^. (Note 
that the ordinate scaling for H2 should be doubled when comparing the two curves in mass.) The 
abscissa is the height from the galactic plane. 

Figure 3. (a) The {R, z) distribution of the HI gas density. Contours are drawn at intervals of 0.01 H 
cm"''. (6) The {R,z) distribution of the H2 gas density. Contours are drawn at intervals of 0.1 H2 
cm""'. 

Figure 4. The {R, z) distribution of the H2 gas density for the whole range of the CO data. Contours 
are drawn at intervals of 0.2 H2 cm~^. 

Figure 5. (a)Radial distributions of HI gas density (doted line) and H2 gas density (solid line) averaged 
in \z\ < 100 pc. (&) The same, but including the radial variation of the CO-to-H2 conversion factor 
C2 as given by equation (5). 

Figure 6. The molecular fraction at I = 53° as a function of z. The ordinate is /moi- 

Figure 7. The {R, z) distribution of the molecular fraction. Contours are drawn at an interval of 0.2. 
Dashed lines show /moi= 0.2, 0.4, and solid lines /moi= 0.6, 0.8. 

Figure 8. (a) Radial variations of the molecular fraction averaged within \z\ < 50 pc (solid line) and 
200 pc (dashed line). (6) The same as (a) but including the effect of radial dependency of CO-to-H2 
conversion factor. 

Figure 9. The molecular fraction calculated theoretically for I = 53°. The solid line is for the case 
where the observed ntot(-z) was used to calculate P{z). The dashed line corresponds to a gaussian 
distribution of P{z). The doted line shows the observed /moi (Fig- 6). 
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Table 1. Observational parameters 





HT 

±±± 


rn (7 — 1 ^ c\] 




Rania and Tjorkman f1Q84l 


KnaoD et al f 19851 


xeiescope 




jjen j-fcLuoiaxoiies 


Aperture 


QH/I TVl 

oU'i m 


7 m 


Wavelength 


z± cm 


2.6 mm 


HPBW 


4' 


1'40" 


TimiP"itnHp P'X'tpnt 


33° 64° 


10° ~ 77° 


Latitude extent 


±3° 


~ ±2° 


Av 


l.OSkms-i 


0.65 km s-i 


Ab 


2' 


2' 


R (kpc)* 


4.68 ~ 7.65 


1.49 ~ 8.29 


Longitude interval 


AsinZ = 0.025, or Al 


= l°.43/cos/ 


Radius interval 


AR = RoAsml = 


212.5 pc 



The distance to the Galactic Center is taken to be Ro = 8.5 kpc. 



11 



